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Abstract 

The completed IceCube Observatory, the first km^ neutrino telescope, is already providing the most stringent limits 
on the flux of high energy cosmic neutrinos from point-like and diffuse galactic and extra-galactic sources. The non- 
detection of extra-terrestrial neutrinos has important consequences on the origin of the cosmic rays. Here the current 
status of astrophysical neutrino searches, and of the observation of a persistent cosmic ray anisotropy above 100 TeV, 
are reviewed. 
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1 Introduction 

One hundred years after their discovery, the origin 
of the cosmic rays is still a mystery. The current 
leading model is that cosmic rays are accelerated in 
diffusive shocks. In this case Supernova Remnants 
(SNRs) in our Galaxy could be the major source of 
cosmic rays up to about 10^^ - 10^'^ eV. The SNR 
energy output in the Galaxy can provide the energy 
budget necessary to maintain the presently observed 
population of galactic cosmic-rays. In particular, in 
order to achieve such high energies it is expected that 
acceleration occurs during the relatively short period 
in the SNR evolution between the end of free expan- 
sion and the beginning of the so-called Sedov phase. 
This period is about 10'^ years from the explosion 
when the shock velocity is high enough to allow for 
efficient acceleration. At energies in excess of about 
10^'^ eV, Active Galactic Nuclei (AGN) and Gamma 
Ray Bursts (GRB) could play an important role in 
the origin of the extra-galactic cosmic rays. 

Since cosmic rays are deflected by magnetic fields, 
it is not possible to associate them to their sources. 
However, if hadronic particles are accelerated, a frac- 
tion of them would interact within their sources or 
in surrounding molecular clouds to produce mesons. 
The mesons eventually decay into high energy 7 rays 
and neutrinos with an energy spectrum ^ E^'^ of 
the accelerated cosmic rays. The remaining hadronic 
particles propagate until their detection on Earth. 
Detection of 7 rays and neutrinos from individual 
galactic or extra-galactic source candidates of cosmic 
rays, or from extended molecular clouds, is there- 

^most probably SNR older than several thousand years no 



fore a method to indirectly probe the origin of cosmic 
rays. 

During the last decade, detection of 7 rays from 
galactic sources has been successfully achieved by 
satellite experiments such as AGILE and Fermi up 
to 10 and 100 GeV, respectively. Imaging Cherenkov 
Telescope Arrays such as MAGIC, VERITAS and 
H.E.S.S., and water Cherenkov detectors such as Mi- 
lagro have made measurements up to 0(10 TeV). 
High energy direct emission from old SNRs appears 
to be inconsistent with hadronic acceleratiorl^. It is 
interesting, however, that delayed secondary 7 ray 
emissions can be produced by the most energetic par- 
ticles that escaped the acceleration region when they 
propagate through molecular clouds that surround 
the star forming regions [1]. With this mechanism, 
indirect evidence of hadronic acceleration is present 
even when SNR are several 10** years old. In fact, the 
detection of an extended emission of TeV 7 rays from 
the Galactic Center by H.E.S.S., which is attributed 
to cosmic rays accelerated by SNR G0.9-f0.1 inter- 
acting with the surrounding clouds, might provide 
the first evidence of hadronic acceleration [5]. The 
most compelling evidence currently comes from low 
energy 7 ray emission from the regions surrounding 
the intermediate-age SNR W44. AGILE observations 
in the energy range of 50 MeV - 10 GeV 3 and Fermi 
observations up to 100 GeV [1] show that while lep- 
tonic models fail to describe simultaneously 7 and 
radio emissions without requiring too large circum- 
stellar densities, the hadronic models are consistent 
with experimental constraints from radio, optical, X 
and 7 rays observations. Although the 7 ray energy 

ger efficiently accelerate cosmic rays. 
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spectrum is consistent with a proton spectral index of 
3 and a low energy cut-off of approximately 10 GeVH, 
the hadronic origin of the observed emission is con- 
sidered likely. The observed steep spectrum and low 
energy cut-off may be caused by suppression of efh- 
cient particle acceleration in the dense environment 
of this source 151 . Ion- neutral collisions in the weakly 
ionized dense gas surrounding the remnant lead to a 
softer spectrum as well as to damping of the plasma 
Alfven waves that form the shock. The resulting poor 
particle confinement leads to a low energy cutoff [6 . 

Other than the specific properties of single ob- 
jects, evidence of an instance of hadronic accelera- 
tion is a very important step towards the discovery 
of the origin of cosmic rays. However, this would 
not mean that all galactic cosmic rays are necessarily 
accelerated in SNR. If cosmic ray acceleration occurs 
predominantly on a larger scale, such as in superbub- 
bles [7] or in the Galaxy cluster medium where par- 
ticles could be accelerated to ultra- high energies 0, 
the search for the origin of cosmic rays should con- 
centrate on extended sources or diffuse fluxes. 

While the TeV 7 ray horizon is limited within our 
Galaxy, because of absorption in the infrared and mi- 
crowave cosmic background, the GeV 7 emissions can 
be observed within about 100 Mpc making it possible 
to search for extragalactic sources of cosmic rays. On 
the other hand, detection of neutrinos from individ- 
ual sources are an efficient and unambiguous probe 
for the high energy hadronic acceleration mechanism, 
and therefore for the sources of cosmic rays. How- 
ever, the very same property that makes neutrinos 
an excellent cosmic messenger also makes them dif- 
ficult to detect. Thus large instrumented volume of 
target matter is required to capture sufficient event 
statistics. 

The IceCube Neutrino Observatory (see Fig. [T]), 
completed in December 2010, is currently the only 
km^ scale neutrino telescope collecting data. The 
observatory consists of an array of 5,160 optical sen- 
sors arranged along 86 cables (or strings) between 
1,450 and 2,450 meters below the geographic South 
Pole, where the antarctic ice is particularly transpar- 
ent. IceCube includes a surface shower array, IceTop, 
and a dense instrumented core with a lower energy 
threshold, DeepCore. The surface array, IceTop, is 
81 stations each consisting of two tanks of frozen 
clean water with each tank containing two optical 
sensors. IceTop, using events in coincidence with the 
deep IceCube array, provides the measurement of the 
spectrum and mass composition of cosmic rays at the 
knee and up to about 10^^ eV. The DeepCore sub- 
array, consisting of 6 densely instrumented strings 
located at the bottom-center of IceCube, lowers the 
observatory neutrino energy threshold to about 10 
GeV. DeepCore uses the surrounding IceCube instru- 

^this is the reason why such a source was not observed at 



mented volume as a veto for the background of cosmic 
ray induced through-going muon bundles, thus en- 
hancing the detection of down-going neutrinos within 
the Deep Core volume. Veto rejection power in ex- 
cess of 10* has been achieved The basic detection 
component of IceCube is the Digital Optical Module 
(DOM) which consists of a 10-inch Hamamatsu pho- 
tomultiplier tube (PMT) and its own data acquisition 
(DAQ) circuitry enclosed in a pressure-resistant glass 
sphere. The DOMs detect, digitize and timestamp 
the signals from Cherenkov radiation photons. Their 
main DAQ board is connected to the central DAQ 
in the IceCube Laboratory at the surface, where the 
global trigger is determined |11| . The construction of 
IceCube started in 2004 and physics quality data tak- 
ing commenced in 2006. With this early data the ob- 
servatory is providing the most stringent limits on the 
flux of high energy neutrinos from extra-terrestrial 
origin, and therefore strong constraints on the models 
of individual sources of cosmic rays and unidentified 
diffuse sources. At the same time, IceCube has accu- 
mulated a large number of cosmic ray induced neutri- 
nos produced in the atmosphere, making it possible 
to probe the combined effect of hadronic interaction 
models, cosmic ray spectrum and composition on the 
neutrino spectrum up to a few hundred TeV [TO] . 




Figure 1: A schematic view of the IceCube Observa- 
tory with the surface array IceTop and the densely 
instrumented DeepCore. 

In the search for high energy neutrinos, the large 
exposure of IceCube makes it possible to collect an 
unprecedented number of events in the form of bun- 
dles of high energy muons generated in the cosmic 
ray induced extensive air showers. Although these 
events represent an overwhelming background in the 
neutrino searches, they make it possible, for the first 
time, to determine the degree of anisotropy of cosmic 
rays from a few TeV to several PeV of particle energy. 
The persistence of a cosmic ray anisotropy at high 

energy. 
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energy raises the question of the responsible mecha- 
nism. The notion that cosmic ray anisotropy might 
be connected to the distribution of nearby and recent 
supernovae is intriguing, and might thus provide a 
new probe into the origin of the cosmic rays. On the 
other hand the complex energy- dependent topology 
suggests that non-diffusive processes in the local in- 
terstellar medium most probably play an important 
role. 

2 Physics Results 

If the signals from detected Cherenkov photons sat- 
isfy specific trigger conditions, an event is defined and 
recorded by the surface data acquisition system. On- 
line data filtering at the South Pole reduces the event 
volume to about 10% of the trigger rate, based on a 
series of reconstruction and filter algorithms aimed 
to select events based on directionality, topology and 
energy [15]. The filter makes it possible to trans- 
fer data via satellite from the experimental site for 
prompt physics analyses. 

2.1 Atmospheric neutrinos 

Of the events that trigger IceCube, the vast majority 
are muon bundles produced by the impact of pri- 
mary cosmic rays in the atmosphere. Only a small 
fraction of the detected events (~10~^) are muons 
produced by the charged current interaction of at- 
mospheric muon neutrinos. The easiest way to re- 
ject the down-going muon bundle background is to 
exclusively select well reconstructed up-going events, 
since these can only be produced by neutrinos cross- 
ing the Earth and interacting in the matter surround- 
ing the detector. Depending on the detector configu- 
ration and on the specific reconstruction algorithms 
and event selection utilized, the atmospheric neutrino 
sample is characterized by a directional resolution of 
better than 1° above 1 TeV. The corresponding reso- 
lution in the estimation of the muon energy is about 
0.2-0.3 (in loglO of the energy) for crossing track- 
like events, and about 0.1 or better for contained 
cascade-like events. Typically, 30%-40% of the up- 
going events survive the selection with a background 
contamination of less than about 1% (see Tab.[T]). 

Table 1 : Mean rate of muon bundles and atmospheric 
neutrinos after final event selection for different string 
configurations of the IceCube Observatory (numbers 
in italic are predictions). 



strings 


year 


moan /i rate 


final i/p rate 


22 


2007 


500 Hz 


18 / day 


40 


2008 


1100 Hz 


40 / day 


59 


2009 


1700 Hz 


130 / day 


79 


2010 


2000 Hz 


no / day 


86 


2011 


2100 Hz 


200 / day 



The atmospheric neutrino sample collected by 
IceCube over the years is the largest ever recorded 



and currently reaches energies near 400 TeV (see 
Fig. [5]). For the first time the precision of this mea- 
surement is providing a powerful tool to constrain 
the effects of high energy hadronic interaction models 
that represent our present knowledge of the cosmic 
ray induced extensive air showers and the spectrum 
and composition of primary cosmic rays |10j . 
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Figure 2: Collection of theoretical calculations and 
experimental measurements of the atmospheric neu- 
trino spectrum. Shown is the predicted conventional 
v^ + v^ (blue line) and Ve + i>e (red line) fiux from 16 , 
and the predicted prompt flux of neutrinos (magenta 
band) from [TTJ. The unfolded energy spectrum [18] 
(black filled circles) and forward folded spectrum [19] 
(gray band) from the 40-string IceCube configura- 
tion, unfolded spectrum [2D] (black open circles) 
and forward folded spectrum [21] (ecru band) from 
AMANDA are presented. The results from Super- 
K [2^ (aqua band) and that from Frejus [23* (black 
filled squares for -\- and black open squares for 
Ve + Ve) are also presented. 

2.2 Search for astrophysical z/'s 

Atmospheric neutrinos represent an irreducible back- 
ground for the search of high energy astrophysical 
neutrinos. If hadronic acceleration is the underly- 
ing process of high energy cosmic ray production 
and 7 ray observations in galactic and extra-galactic 
sources, the charged mesons could produce enough 
neutrinos to be observed in a detector the size of Ice- 
Cube. Fig. [3] shows the sensitivity (90% CL) of Ice- 
Cube for the full-sky search of steady point sources 
of E^^ muon neutrinos as a function of declination, 
along with that of other experiments. The extension 
of the point source search to the southern hemisphere 
is made possible by a high energy event selection that 
rejects the background down-going events by five or- 
ders of magnitude, and restricts neutrino energies to 
above 100 TeV. Still dominated by high energy large 
muon bundles, this makes the southern hemisphere 
poor in atmospheric neutrinos yielding a low neutrino 
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detection sensitivity. Nevertheless, this provides Ice- 
Cube with a full-sky view that complements cover- 
age of the neutrino telescopes in the Mediterranean. 
The figure shows the sensitivities from IceCube and 
other observatories (interpreted as the median upper 
limit we expect to observe from individual sources 
across the sky) along with upper limits from selected 
sources. The sensitivity is reaching the level of cur- 
rent predictions for flux from astrophysical sources 
{i.e. below lO^^x E'^ TeV cm'^ s"^) although the 
discovery potential, defined to be 5a for 50% of the 
trials, is typically a factor of three higher than the 
sensitivity. Therefore constraints on the parameters 
of hadronic acceleration models are starting to de- 
velop. 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 
sin{5) [degrees] 



Figure 3: Sensitivity (90% CL) for a full-sky search 
of steady point sources of muon neutrinos with an 
energy spectrum as a function of declination an- 
gle for IceCube and other experiments. Note that for 
IceCube, events with 5 < 0° are down-going, coming 
from the southern hemisphere, and events with 5 >0° 
are up-going and come from the northern hemisphere. 
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Figure 4: Upper limits (90% CL) for neutrino 
searches in coincidence with Gamma Ray Bursts with 
40 strings of IceCube and the combined 40- and 59- 
string detector configurations |32j . Also shown is the 
Waxman & Bahcall predicted average flux j33j. 

Searches for neutrinos from transient |30| and pe- 
riodic [21] sources have also been performed. In 
particular, a time window scan for transient sources 
(with no external triggers) shows that the discovery 
potential drops by a factor of 2 if searching for 1 



day duration flares. A particular search for transient 
sources is that for neutrinos from GRB. For the first 
time, the IceCube Observatory has provided a defini- 
tive test of the GRB models with the most stringent 
constraints. Fig. |4] shows the upper limits obtained 
with the data collected by the 40-string configura- 
tion of IceCube and by the combined data of the 
40- and 59-string configurations [32]. For each de- 
tector configuration, a list of GRBs detected during 
the corresponding physics runs was compiled and the 
predicted neutrino flux was calculated based on the 
7 ray spectrum shown in |34j . The corresponding 
stacked neutrino flux was used to search for events 
collected within the time window in which 5% to 
95% of the fluence is recorded {i.e. Tgo)- The upper 
limit is about 3 times below the predicted fiux of the 
Waxman & Bahcall model, challenging the hypoth- 
esis that GRB are the sources of Ultra High Energy 
Cosmic Rays (UHECR). This result has profound 
consequences for the predicted flux of neutrinos pro- 
duced by the interaction of UHECR with the cosmic 
microwave background, the so-called cosmogenic neu- 
trinos, as well as for the GeV- TeV 7 ray background 
flux (see for instance [3Sl[3n])- It is important to note 
that it was recently shown that the fireball model 
with refined assumptions yields a 10 times smaller 
predicted flux (see [37l|38]). There is the possibility 
that the bulk of cosmic rays does not originate from 
individual sources, but from large-scale acceleration 
processes in superbubbles or even Galaxy clusters. 
In addition, unresolved sources of cosmic rays over 
cosmological times are expected to have produced 
detectable fluxes of diffuse neutrinos. Since shock 
acceleration is expected to provide an ~E~^ energy 
spectrum, harder than the ~E~'^-^ of the atmospheric 
neutrinos, the diffuse fiux is expected to dominate 
at high energy where the sensitivity is strongly de- 
pendent on the experimental quality of the selected 
events. Fig. [S] shows a collection of sensitivities and 
upper limits (90% CL) for an E~^ flux of i/^ + D^i, 
from AMANDA, Antares and various IceCube con- 
figurations compared to the experimental and theo- 
retical flux of the atmospheric neutrinos and various 
models of astrophysical neutrinos. The most recent 
results lie below the Waxmann & Bahcall neutrino 
bound [32, again indicating IceCube's potential for 
discovering the origin of cosmic rays. 

In the Ultra High Energy range (UHE), above 
~10^ GeV, IceCube is reaching a competitive sen- 
sitivity as well. At this level one begins to reach 
current models of cosmogenic neutrino production 
(see Fig. [6]) that are simultaneously constrained by 
the current observations of UHECRs and the GeV 7 
rays by Fermi-LAT [49]. Taking into account that 
UHECR mass composition is a key ingredient for 
the absolute flux and spectral shape of cosmogenic 
neutrinos [35], its large uncertainty still weighs pro- 
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foundly on cur.-ent models. This means that although 
the IceCube sensitivity to UHE neutrinos is currently 
the best ever achieved below 10^" GeV it might be 
still far from the actual flux. From this point of 
view, the current developments toward a radio ar- 
ray in Antarctica, such as Askaryan Radio Array 
(ARA) [5T] is a natural extension toward the highest 
energies. 



on a background of conventional atmospheric neutri- 
nos of 0.3. The events deposited an energy in the de- 
tector of about 1 PeV, and further study is underway 
to determine their nature. One possible hypothesis 
is that these events represent an upper fluctuation of 
the prompt neutrino production in the atmosphere 
from the decay of heavy charm mesons. 
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Figure 5: Experimental upper limits (90% CL) for 
the diffuse muon neutrino flux (including the prelim- 
inary result from the 59-string configuration of Ice- 
Cube) along with atmospheric neutrino observations 
and theoretical models of atmospheric and extra- 
terrestrial neutrino fluxes. From top to bottom in 
the legend [32 SQl SIl [111 [161 Eg |43l Ell Sg 
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Figure 6: Preliminary sensitivity (90% CL) for the 
detection of UHE neutrinos, compared to other ex- 
perimental results and to predictions [48l|49l[5^. The 
sensitivity curves are evaluated at each decade of en- 
ergy 

It is worth noting that the preliminary sensitiv- 
ity for an arbitrary spectrum, shown in Fig. [51 has 
a minimum just above 1 PeV, where no significant 
cosmogenic neutrino flux is expected. In the exper- 
imental analysis performed on data collected during 
2010-12, where events with a large number of de- 
tected photons were selected, two events were found 



2.3 Cosmic ray anisotropy 

The large number of muon bundle events collected 
by IceCube (about lO^'^-lO^^ each year, depending on 
the detector configuration) makes it possible to study 
the arrival direction distribution of the cosmic rays 
at a level of about 10~^. The bundles of highly col- 
limated atmospheric muons share the same direction 
as the parent cosmic ray particle. Since this study 
does not require highly well reconstructed muon di- 
rections, all collected and reconstructed events with 
a median angular resolution of about 3° are used. 
Using full simulation of cosmic ray induced extensive 
air shower we find that the median particle energy of 
the IceCube data sample is about 20 TeV. With these 
data IceCube provides the first high statistics deter- 
mination of the anisotropy of galactic cosmic rays 
in the southern hemisphere in the multi-TeV energy 
range. 

The large scale anisotropy observed by Ice- 
Cube [52 1 appears to complement the observations 
in the northern hemisphere, providing for the first 
time an all-sky view of TeV cosmic ray arrival direc- 
tions. The sky map obtained by subtracting an aver- 
aged map (over a scale of 30°-60°) from the data [53], 
shows significant small angular scale structures in the 
cosmic ray anisotropy, similarly to observations in 
the northern hemisphere |54[ 155) . Another interest- 
ing result obtained by IceCube is the persistence of 
the anisotropy at an energy in excess of 100 TeV. At 
such energies a different structure is observed that 
can be interpreted in terms of a different phase |56] 
as already reported by the EAS-TOP shower array in 
the northern hemisphere for the first time [58]. The 
observation at high energy was recently confirmed 
by the preliminary result from the IceTop shower 
array [57]. The change of the anisotropy pattern 
at about 100 TeV may suggest that the heliosphere 
could have an effect in hipping the apparent direc- 
tion of the anisotropy. In fact, at about 100 TeV 
the cosmic rays' gyro-radius in the 3 /iG local in- 
terstellar magnetic field is of the order of magnitude 
of the elongated heliosphere. Below this energy scale 
the scattering processes on the heliospheric perturba- 
tions at the boundary with the interstellar magnetic 
field might be the dominant processes affecting the 
global cosmic ray arrival distribution and the small 
angular structure as well (see [59] where a review of 
other proposed models is also given). The Milagro 
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observation of a likely harder than average cosmic 
ray spectrum from the localized excess region toward 
the direction of the heliotail, the so-called region B 
in [54] and also observed by ARGO-YBJ shower ar- 
ray [55) . have triggered astrophysics interpretations 
(see [601 EU |62]). However, this may suggest that 
some type of re- acceleration mechanism associated 
with cosmic ray propagation in the turbulent helio- 
spheric tail might occur [63l|64]. On the other hand, 
the TeV cosmic ray anisotropy is a tracer of the lo- 
cal interstellar magnetic field, and it might indicate 
cosmic ray streaming along the magnetic field lines 
due to the Loop I shell expanding from the Scorpion- 
Centaurus Association |65j . 

If the local propagation effects on the cosmic ray 
anisotropy below 100 TeV are dominant, at higher 
energy it is reasonable to believe that the persistent 
anisotropy might be a natural consequence of the 
stochastic nature of cosmic ray galactic sources, in 
particular nearby and recent SNRs [66l [67l [68] . 
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